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We present cosmological constraints on deviations of Newton’s constant at large scales, analyzing 
latest cosmic microwave background (CMB) anisotropies and primordial abundances of light ele¬ 
ments synthesized by big bang nucleosynthesis (BBN). BBN limits the possible deviation at typical 
scales of BBN epoch, say at 10® ^ lO’^^ m, to he between —5% and -|-1% of the experimental value, 
and CMB restricts the deviation at larger scales 10^ ~ 10® pc to be between —26% and -1-66% at 
the 2(7 conhdence level. The cosmological constraints are compared with the astronomical one from 
the evolution of isochrone of globular clusters. 


PACS numbers: 98.80.Cq, 98.65.Dx, 98.70.Vc 

Newton’s law of gravitation has been extensively tested 
and verified in three length scales: the laboratory scales 
r < 1 m Q, the geophysical scales r k, 100 m Q, and 
the astronomical scales r « 10® m 1^. Such measure¬ 
ments nicely agree with the inverse square law within 
their experimental or observational uncertainty ill In 
particular, the first two measurements at the laboratory 
and geophysical scales succeeded also in determining the 
experimental value Gn of Newton’s constant, and the 
value determined at such terrestrial scales is applied for 
all phenomena from Planck scale to cosmological scale. 

The astronomical measurements Q, mainly through 
planetary and satellite orbits, yield a strong constraint 
on the deviation from the inverse square law. However, 
it should be noted that the measurements can not give 
any information about the value of Newton’s constant G 
itself without evaluating masses M of interacting bodies, 
since constraint is possible only on GM. Therefore, the 
measurements can not exclude the possibility of differ¬ 
ent value of G at astronomical and cosmological scales, 
if G is almost constant at limited scales relevant to the 
measurements. In particular, we have only a poor knowl¬ 
edge at scales larger than the solar system, say r > 1 pc 
R:! 3 X 10^® m. s Interesting trials to solve this problem 
were recently reported HQ , in which the deviation of G 
at Mpc scales is restricted by the power of the clustering 
of galaxies . 

The possibility that Newton’s constant at laboratory 
scale, Gn, is different from that at very large scales. Goo, 
arises in many context. Historically, studies toward the 
problem of unifying gravity with the other fundamental 
forces suggested a departure from Newtonian gravity in 
the range 10 — 100 m. It is often assumed that 

such a correction can be represented by the addition of 
Yukawa term to the conventional gravitational potential: 
V = for G(r) = Goo (l + , where a is 

the relative weight of the non-Newtonian term. In this 
expression, at cosmological distances r satisfying r ^ A, 


the exponential term vanishes, so that G(r) = Goo- On 
the other hand, for r of experimental scales which satisfies 
r <C A, the exponential becomes unity and G(r) recovers 
Gn, that is, Gn = Goo (1 -I- a) . 

Recently several types of higher-dimensional theories 
of gravity, motivated by superstring, have been proposed 
and many researchers pay great attention to the extra di¬ 
mension scenario. As a characteristic feature, all the the¬ 
ories lead to deviations from the conventional Newton’s 
law 11113, since the theories allow graviton to propa¬ 
gate in higher-dimensional spacetime. Among them, an 
interesting idea was proposed by Dvali, Gabadadze and 
Porrati. In the model the present accelerating expan¬ 
sion of the universe is attributed to leaking gravity into 
extra dimension El This idea reproduces the present 
cosmic acceleration without dark energy component, and 
consequently predicts the modification of Newton’s law 
at cosmological scales. Another interesting proposal is a 
braneworld model with Gauss-Bonnet term, which sug¬ 
gests Gn = q^Goo with a model parameter a [T^ . 

These theoretical suggestions indicate that it is quite 
important to place possible constraints on the value of G 
at astronomical and cosmological scales. In this paper, 
we take a simple parameterization of 

G(r) = CGn , (1) 

for a finite range of r relevant to the measurement which 
we are considering. As such a measurement, we consider 
CMB anisotropies and primordial abundances created at 
BBN epoch, and put the constraints on the value of ^ at 
two different scales relevant to these observations. The 
cosmological constraints are compared with an astronom¬ 
ical constraint determined from the isochrone of globular 
clusters. 

The observed primordial light-element abundances 
constrain the value of G during the BBN epoch from the 
time of weak reaction freezeout (t ^ 1 sec, T ~ 1 MeV) 
to the freezeout of nuclear reactions (t 10^ sec, T ~ 
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10 keV). In this epoch, the length of cosmic horizon 
varies from 10® m to 10^^ m, and thus BBN can con¬ 
strain Newton’s constant at these scales. 

The primordial helium abundance is obtained by mea¬ 
suring extra galactic HII regions. We adopt range of 
Yp = 0.2452 ± 0.0015 for the helium abundance. 
The primordial deuterium is best determined from its 
absorption lines in high redshift Lyman a clouds along 
the lines of sight to background quasars. For deuterium 
there is a similar possibility for either a high or low value. 
For the present discussion, however, we shall adopt the 
generally accepted low value for the D/H abundance, 

D/H = 2.78t!^;^^ x 10”® Hi- 

The increase of Newton’s constant causes the increase 
of the universal expansion rate. This makes the neutron- 
to-proton ratio larger, because the weak reactions freeze- 
out at a higher temperature, and also because there is 
less time for neutrons to decay between the time of weak- 
reaction freezeout and the onset of BBN. Consequently, 
larger value of G during BBN epoch yields a larger "^He 
abundance, since most of the free neutrons are converted 
into “^He nuclei. D/H also increases largely because the 
reactions destroying deuterium fall out of nuclear statisti¬ 
cal equilibrium while the deuterium abundance is higher. 
Similarly, there is less time for the destructive reaction 
^Li(p, a)'^He. This causes ^Li to be more abundant for 
ry < 3 X 10“^°. However, there is also less time for the 
^He(®He, 7 )^Be reaction to occur. This causes ^Li to be 
less abundant for ry > 3 x 10“^° Ea 

The upper limit of Newton’s constant comes from ^He 
upper bound and D/H upper bound. The lower limit 
comes from the lower bounds. We note that the con¬ 
straint from ^Li is not consistent with those from ^He 
and D/H, even when we vary In the present analy¬ 
sis, however, we omit the constraint from ^Li abundance, 
since it involves an uncertainty more largely than the 
other primordial elements do. The BBN constraint thus 
obtained is 0.95 < C < 1.01. 

The temperature fluctuations at recombination ob¬ 
served through the CMB anisotropies contain much in¬ 
formation on many kinds of cosmological parameters and 
evolution of perturbations at a wide range of scales. Typ¬ 
ically, the scale which can be explored by CMB obser¬ 
vations currently available is from the horizon scale at 
present (~ Gpc) to ~ 10 Mpc in comoving coordinate. 
This shows that the scales relevant to CMB are from 10^ 
pc to ~ Gpc, if we consider the evolution of perturbations 
from horizon crossing of each Fourier mode; for example, 
10^ pc is the horizon scale at the time when the mode 
of ^ 10 Mpc in comoving coordinate enters the horizon. 
Thus it follows that CMB can constrain the value of ^ at 
scales larger than ~ 10^ pc. Here, in order to calculate 
CMB anisotropies in a consistent manner we assume that 
the scale dependence of Newton’s constant is very weak 
at the relevant scales, which is consistent with a simple 
parameterization of Eq. (1). 

In order to obtain a constraint on ^ from latest CMB 
anisotropy data sets, we generate CMB angular power 
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FIG. 1: Predicted BBN light-element abundances v.s. the 
baryon-to-photon ratio rjio- “^He, D and ^Li abundances are 
shown in the top, center and bottom panels, respectively. 
They are compared with the observationally inferred primor¬ 
dial abundances (horizontal lines). Plotted are models with 
5 = 1.01 (solid) and ^ = 0.95 (dashed). In our analysis neu¬ 
tron life time is taken to be the average value of t„ = 878.5 
sec and Tn = 885.7 sec M- 






FIG. 2: Gonstraints from the primordial abundances (lines) 
and CMB (contours) on ryio and C The shaded region de¬ 
notes parameters allowed by BBN. The contours show the 
marginalized 1, 2cr limits on this parameter plane from hts to 
the CMB power spectrum. The solid and dashed contours 
correspond to the limits from WMAP data alone, and those 
from WMAP, CBI and ACBAR data sets, respectively. 


spectra Ci in a wide range of ^ by using a Boltzmann code 
of CMB FAST It is well known, however, that in ad¬ 

dition to ^ there exist many other cosmological param¬ 
eters relevant to CMB. Thus, we explore the likelihood 
in seven dimensional parameter space, i.e., (baryon 

density), (cold dark matter density), h (Hubble pa¬ 
rameter), Zre (reionization redshift), Ug (power spectrum 
index), Ag (overall amplitude), and We then marginal¬ 
ize over nuisance parameters through the use of Markov 
Chain Monte Carlo technique [T^ . 
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FIG. 3: CMB power spectrum with and without the variation 
of Higher peaks are more severely damped as ^ increases, 
while the height of the hrst peak is almost unchanged. 


FIG. 4: Marginalized probability distribution of The solid 
and dashed lines correspond to the probability distributions 
obtained by WMAP data alone and by WMAP, CBI and 
AGBAR data sets, respectively. 


The most distinguishable effects of changing Newton’s 
constant appear at the amplitude of the acoustic peaks 
in the CMB power spectrum as shown in Fig. 0 The 
main reason for this is that, as already found in | 2 (l| . 
the visibility function, g^r) = kexp(—/c), changes with 
where k is the optical depth of the Thomson scat¬ 
tering. More specihcally, increasing Newton’s constant 
makes the expansion of the universe faster at a give red- 
shift, and makes it more difficult for proton and electron 
to recombine to form hydrogen atom. This leads to larger 
ionization fraction and broader visibility function at last 
scattering epochs, which damp the anisotropies at small 
scales due to canceling effect. 

Second, in addition to the effect discussed above, we 
find that increasing Newton’s constant suppresses the 
second and higher acoustic peaks even larger since for 
the increase of ^ the scale of the first acoustic peak 
(~ tdec OC is shifted more largely than the dif¬ 

fusion scale for photons to spread through the random 
walk does, (~ oc Thus, the shape of acoustic 

peaks can be used to constrain the variation of 

Figure 01 shows the marginalized likelihood of We 
obtain from the figure that 0.74 < ^ < 1.66 by WMAP 
data alone [HI, 0.75 ^ < 1.74 by WMAP, CBI and 
ACBAR data sets [^ [^ , at 95% confidence level. 

Another constraint on the value of Newton’s constant 
can be obtained by analyzing the age of stars in globular 
clusters. The key idea is that increasing Newton’s con¬ 
stant causes stars to burn faster Thus, this allows 
us to constrain ^ at stellar scale ~ 10 ® m, as we shall see 
below, by analyzing the timing of main sequence turn off. 

Let us assume that luminosity of star depends on 
Newton’s constant G and helium abundance T, approx¬ 
imately as L oc y{Y)g{G), where y and g are functions 
of Y and G [^ . Since helium production should be pro¬ 
portional to the luminosity, we have ^ oc y{Y)g{G). A 
star which departs from the main sequence today (to) 
should be considered to have Y Ki 1 a.t its center so 


that Jy. oc g{G) ^ dt. We further assume that 

g{G) oc G'*', where ^ » 5.6 have been obtained from nu¬ 
merical simulation [^. From the fact that the l.h.s. of 
the above equation does not depend on G and time, we 
have the relation, 

n = C f dt = . ( 2 ) 

J to—T 

Here r* is the apparent turn-off age, which should be 
obtained by analyzing HR diagram of a globular cluster 
with the standard value of G, and r is the true age of the 
globular cluster. Thus, if information on the true age of 
globular cluster is available, the globular cluster can be 
used to constrain ^; 

( 3 ) 

If we take r^ax = 15.8 (2 a upper bound on the expansion 
age of the universe obtained by our analysis in Sec. HI 
including the variation of ^), and conservatively assume 
that Tniin = 10 Gyr j^, we then obtain 

0.93 < C < 1.09 , (4) 

where we use 7 = 5.6 and t, = 12.9 ± 2.9 Gyr, which is 
age of the galactic globular clusters [^ . 

All the higher-dimensional theories of gravity proposed 
recently allow Newton’s constant to be scale-dependent. 
In this paper, assuming the dependence is weak for hori¬ 
zon scales in BBN epoch (10® ^ 10^® m) and also for 
those in GMB epoch (10® ^ 10® pc), we place constraints 
on ^ = G/Gn at the cosmological scales. An important 
point is that the present analysis yields constraints on 
the value of G itself, while other astronomical tests of 
the inverse square law do so only on the value of GM 
including unknown mass M of interacting bodies. 
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FIG. 5: Limits on the variation off from various observations. 
Limit from SDSS is taken from [g. 
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Increasing Newton’s constant enhances the universal 
expansion rate, and then leads to larger helium and deu¬ 
terium abundances produced at BBN epoch. We have 
re-examined this effect includin g th e latest experimen¬ 
tal data on the neutron lifetime [Ig . We found that 
the experimental value Gn (C = 1 ) is now quite consis¬ 
tent with the observed abundances of primordial light el¬ 
ements, and the variation of Newton’s constant is tightly 
constrained to 0.95 < ^ < 1.01. 

The variation of Newton’s constant also affects the 
power spectrum of CMB anisotropies through the change 
of the recombination and photon diffusion processes. We 
found that the difference emerges at smaller scales. Thus, 


observations at higher multipoles are essential to put a 
tighter constraint. However, even when higher multi¬ 
poles data currently available from CBI and ACBAR are 
included, we found no improvement in constraint on f, 
because of scatters in data at higher multipoles. WMAP 
data alone place a constraint: 0.74 < ^ < 1.66. If we 
combine CBI and ACBAR data sets, the constraint be¬ 
comes 0.75 < ^ < 1.74. 

In Fig. 13 we summarize results of the current work. 
The value of ^ is fixed to one at laboratory scale ~ I m 
by direct experiments. We have two possibilities of tran¬ 
sition from short distance regime where G = Gn to long 
distance one where G = ^Gn', one is geophysical scale 
(i.e., ~ I km - 100 km, where the constraints on the in¬ 
verse square law are relatively weak), the other is scale 
beyond the solar system (> 10^^ m), where we have 
only poor knowledge on G. If we consider the former 
case, BBN gives the tightest constraint on Globu¬ 
lar cluster also gives a consistent but weaker constraint. 
On the other hand, if we consider the latter case, CMB 
anisotropies and galaxy clustering HQ are the only ob¬ 
servations to put constraints on Thus, higher precision 
measurements of CMB anisotropies, particularly in its 
higher multipoles are highly expected to determine the 
value of G at large scales beyond the solar system and 
then to confirm the necessity of the higher-dimensional 
theories of gravity. 
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